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Abstract 

The Weyl-Wigner-Groenewold-Moyal formalism of deformation quantization is applied to cos- 
mological models in the minisuperspace. The quantization procedure is performed explicitly for 
quantum cosmology in a flat minisuperspace. The de Sitter cosmological model is worked out 
in detail and the computation of the Wigner functions for the Hartle-Hawking, Vilenkin and 
Linde wave functions are done numerically. The Wigner function is analytically calculated for 
the Kantowski-Sachs model in (non)commutative quantum cosmology and for string cosmology 
with dilaton exponential potential. Finally, baby universes solutions are described in this context 
and the Wigner function is obtained. 
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I. INTRODUCTION 



One of the fundamental problems of cosmology is the initial singularity. It is commonly 
believed that near this singularity the physical evolution is governed by quantum mechanics. 
In the quantum cosmology framework the whole universe is represented through a wave 
function satisfying the Wheeler-DeWitt equation pLj. It is well known that it represents 
the low energy approximation of string theory, however it contains some non-trivial leading 
order information. The development of quantum cosmology started at the beginning of the 
80's of last century and one of its main ideas is that the universe could be spontaneously 
nucleate out from nothing [2"HTU]. After nucleation the universe can enter into a phase of 
inflationary expansion and continues its evolution to the present time. However, there are 
several important questions that remain to be solved like the general definition of probability, 
time and boundary conditions [IT] . In order to find a unique solution of the Wheeler- 
DeWitt equation it is necessary to impose boundary conditions. In the case of quantum 
mechanics there is an external setup and the boundary conditions can be imposed safely, 
however in 4-dimensional quantum cosmology there is nothing external to the universe and 
the question of which one is the correct prescription for the boundary condition of the 
universe is controversial [12]. There are several proposals for the correct boundary conditions 
in quantum cosmology, for example, the no-boundary proposal of Hartle and Hawking [7] , the 
tunneling proposal of Vilenkin [TU] and the proposal of Linde [E]. In a recent development 
in quantum cosmology a principle of selection in the landscape of string vacua has been 
proposed in [13] , in the context of the minisuperspace. Moreover in [H] different cosmologies 
were defined in terms of a wave function on a compact worldsheet of system of parafermions. 
Much of the information is encoded in the wave function as function of the moduli space. 

Another important issue is how to extract information of the Wheeler-DeWitt equation. 
In general, the configuration space used in quantum cosmology is an infinite dimensional 
space called superspace and it is not amenable to work with. In the study of homogeneous 
universes the infinite dimensional space is truncated to finite degrees of freedom, therefore 
obtaining a particular minisuperspace model. The reduction of superspace to minisuper- 
space is not a rigorous approximation scheme, however there is the hope that minisuperspace 
maintains some of the essential features of quantum cosmology. In classical cosmology ho- 
mogeneity and isotropy are fundamental to describe the universe at large scale and therefore 
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it is expected to have a minisuperspace description in quantum cosmology. 

At very early times of the universe, even before the Planck time, the universe should be 
described by means of quantum cosmology that could take into account effects from string 
theory, supersymmetry and noncommutativity. String cosmology [15], quantum wormholes, 
baby universes and supersymmetric quantum cosmology has become very intensive research 
areas in quantum cosmology. Quantum wormholes [16] are instanton solutions and are 
important in the Euclidean path integral formulation. In the third quantization approach, 
which is an adequate description of topology change in the path integral quantization, the 
wave function is transformed into a quantum field operator which includes operators that 
create and destroy the so called baby universes. 

Supersymmetric quantum cosmology is one of the main research areas [TT1 - I20"] . At the 
time of quantum creation of the universe it is possible that supersymmetry would not yet 
be broken and therefore could have very important consequences in the evolution of the 
universe. Other very interesting effects that could arise in the very early universe are the 
effects of noncommutativity |21j . 

At the end of the late 1980s and early 1990s quantum decoherence, the transition from 
quantum physics to classical physics, was an active research area in quantum cosmology 
[22] [23] (for more recent reviews see, [2U 125]). In Refs. [261430] . quantum cosmology is 
developed in the phase space and the use of Wigner function is showed to be a very useful 
approach to study decoherence. In fact, quantum mechanics in phase space is an appropriate 
formalism to describe quantum mechanical systems (for a review see [SI])- The description 
using the Wigner function has been of considerable interest and usefulness in quantum optics, 
nuclear and particle physics, condensed matter, statistical physics, etc. In particular, the 
description of semiclassical properties and the analysis of the classical limit is more clear in 
the Wigner function formalism. The classical limit from quantum cosmology was studied 
also in [301 E2] , where the use of Wigner function and quantum mechanics in phase space 
is fundamental. It has been argued that quantum decoherence is achieved if the density 
matrix is coarse grained, i.e. averaged over configuration or phase space variables, or in 
an alternative way, if the system interact with some environment which is not monitored. 
However the existence of classical correlations is another characteristic of the classical limit 
and requires the presence of sharp peaks of the Wigner function, but a coarse graining 
produces a spreading of the distribution in phase space. The former arguments demand a 
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subtle coarse graining in which there is a delicate balance between the existence of classical 
correlations and decoherence. A coarse graining could be modeled by a Liouville equation 
with friction and diffusion terms [30J. Some other further results are described in Refs. 

[331 EH- 

The quantum mechanics in phase space is only one part of a complete and consistent 
type of quantization termed: deformation quantization. In this paper we formulate quantum 
cosmology in terms of deformation quantization and rewrite some results of [301 132] in the 
context of this formalism. We will assume that the superspace (of 3-metrics) is flat and 
a Fourier transform can be defined. This is fairly valid for the ample set of examples in 
two dimensions that we present in this paper. The case of curved minisuperspace will be 
discussed in a future communication. 

The deformation quantization formalism is an alternative approach to the canonical and 
path integral quantizations and has its origins in the seminal works by Weyl [35], Wigner 
[3E], Groenewold [37], Moyal [3H], and Vey [32], and is based on the idea of treating with 
quantum mechanics on the phase space. In 1978 Bayen, Flato, Fronsdal, Lichnerowicz 
and Sternheimer [10] introduced its final form in which the quantization is understood as 
a deformation of the classical algebra of observables instead of a change in the nature of 
them. This quantization arises as a deformation of the usual product algebra of the smooth 
functions on the classical phase space and then as a deformation of the Poisson bracket 
algebra. The deformed product is called the *— product which has been proved to exist 
for any symplectic manifold |4~I] . [12] and more recently shown by Kontsevich [13] that it 
also exists for any Poisson manifold. These results in principle allow us to perform the 
quantization of arbitrary Poissonian or symplectic systems and to obtain in an easier way 
the classical limit due to the nature of the *— product. Our aim is to introduce this formalism 
in quantum cosmology and apply the results to simple models. 

The structure of the paper is as follows. In section II we survey the canonical Hamiltonian 
formalism of general relativity and its canonical quantization. In section III we construct the 
Stratonovich-Weyl quantizer, the star-product and the Wigner functional. Section IV is de- 
voted to apply the deformation quantization procedure to several minisuperspace models and 
to obtain their Wigner function. In subsection IV.A we treat de Sitter model and calculate 
numerically the Wigner function for the Hartle-Hawking, Vilenkin and Linde wave functions. 
In subsection IV. B we calculate analytically the Wigner function for the Kantowski-Sachs 
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model in (non) commutative quantum cosmology and we interpret the results. In subsection 
IV. C, we obtain the Wigner function for string cosmology with dilaton exponential potential 
in terms of the Meijer's function. Besides, in subsection IV. D the Wigner function for the 
baby universes solutions is calculated by means of the Moyal product and the annihilation 
and creation operators. Finally, in section V we give our final remarks. 

II. CANONICAL FORMALISM OF GENERAL RELATIVITY 

In this section we briefly overview the Hamiltonian formalism of general relativity. Our 
presentation will not be complete and detailed information can be found in Ref. [TJ EH] . 
Our intention is only to introduce the notation and conventions for future reference along 
the paper. 

A. Hamiltonian Formalism 

We start from the ADM decomposition of general relativity and consider a pseudo Rie- 
mannian manifold (M, g) which is globally hyperbolic. Thus, spacetime M can be decom- 
posed as M = E x R where E is an spatial hypersurface and the metric of a foliation is: 
ds 2 = g^ u dx^dx u = -(N 2 - WN^dt 2 + 2Nidx l dt + hijdx'dx^ with signature (-, +, +, +). 
Here hij is the intrinsic metric on the hypersurface E, N is the lapse function and N l is the 
shift vector. 

The space of all Riemannian 3-metrics and scalar matter configurations $ on E is the so 
called superspace Riem(E) = {hij(x), G E}. Let us denote the space of Riemannian 

metrics on E as Met(E) = {hij(x)\x G E} which is an infinite dimensional manifold. The 
moduli space M., is then defined as the configuration space (superspace) modulo the group 
of diffeomorphisms Dif f (E) i.e. M. = jjjjfpl or for pure gravity M. = • 

The corresponding phase space (Wheeler's phase superspace) V* = T*Met(E) is given 
by the pairs T* = {(hij(x), tt^(x))}, where 7r*- J = 9 q^ h and Leh is the Einstein-Hilbert 
Lagrangian. In the following we will deal with fields at the moment t — (on E) and we 
put hij(x, 0) = hij(x) and Tt^(x,0) = 7r u (x). 

The dynamics of general relativity coupled to matter in this foliation is encoded in the 
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variation of the following action 

1 



S 



dtL 



[ d A x^g~ ( 4 R(g) - 2A) + 2 I d 3 xVhK 

JM JdM 



(1) 



167rGV 

where A R{g) is the scalar curvature in four dimensions depending on the spacetime pseudo- 
Riemannian metric g^, g is the determinant of g^, Gjy is the gravitational constant in N 
dimensions, A is the cosmological constant, K is the trace of the extrinsic curvature in- 
compatible with hij and S m is the matter action. For the case of scalar field matter subject 
to a potential V(Q) we have that 



S m = d A x^g ( - X -g^d^d^ - V($fj . 
The action ([I]) can then be written as 

S I dtd 3 x(-K°N + tt^ - NM± - N ir H 



(2) 



(3) 



where Hj_ and Hi are given below. Thus the canonical Hamiltonian is given by 



d 3 x 7c°N + n % Ni + + 7r$$ 1 — L 



d 3 x 7i° N + 7T l Ni + NHi + N l Ui 



c 



(4) 

where n° = and n l = J^- and ?r$ = |j. The corresponding equations of motion for N 
and Ni yield the Hamiltonian constraint and the momentum constraint, respectively 

-7T 2 - 7T>A + Vk 3 R 



An 2 G ljkl ^n kl - ^( 3 R - 2A) + \^h{^ + ti j $j$j + 2V{$f) = 



Vh 
2^2 



V h 

(G? - 2k 2 T?) = -27rfJ + //' ;< I\/ -.I. = 0, 



(5) 



(6) 



where k 2 = 47tGat, 3 R(h) is the scalar curvature of E, Gijki = \h~ l ^ 2 (hikhji + hiihjk — hijhki)-, 
\j denotes covariant derivative with respect to hij and h is its determinant. 
In this way the Poisson bracket between hij and ir hl is given by 



(7) 



which is one of the most important structures for quantization. 
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B. Canonical Quantization 



The most employed formalism to quantize a physical system is the canonical quantization 
which can be applied in general relativity In this section we describe briefly the general 
procedure to obtain the quantum equations. In the /^-representation, the usual promotion of 
canonical coordinates hij(x), and 7r iJ , 7r$ to the operators can be done in the following 
form: h^hij,®) = h^h^,®}, n lj \hij, $) = -^^rp|^, = ®( x )\ h ij,®) and 

7r<5>\hij, $) = — ihg^n^ \hij, These operators satisfy the commutation relations 



[hij(x),9 



ki 



ih 



M + 6t6l)6{x-v). 



The constraints ^ and ^ have to be imposed at the quantum level in the form 



Ui\*) = 0, 



Ui\V) = . 



(8) 



(9) 



In the coordinate-representation we have 



4k 2 G, 



8 2 



ijkl 



ShjjSh 



+ 



4k 2 V 



^i?0) + 2A + 4^T ! 



^•,$] =0 



(10) 



where (%,$|#) = $] and f 00 = -^^2 + + This constraint is called 



the Wheeler-DeWitt equation. For pure gravity we have 

5 2 Vh 



An 2 Gijki- 



+ 



(- 3 R{h) + 2A) 



0. 



(ii) 



''ShijShki ' 4k 2 

where = ^[/ly]. In the general case the Wheeler-DeWitt equation is not amenable 

to extract physical information of the system. In order to obtain useful information it is a 
common practice to reduce the number of degrees of freedom. In the following sections we 
will consider models with one and two degrees of freedom. 



III. DEFORMATION QUANTIZATION OF WHEELER'S PHASE-SUPERSPACE 

In this section we work in a flat superspace with flat metric Gijki- This case will allow 
us to introduce the Weyl-Wigner-Groenewold-Moyal (WWGM) formalism for gravitational 
systems in a direct way since the existence of the Fourier transform is warranted. The more 
complicated cases of curved (mini)superspaces will be left for a future work. The deformation 
quantization of gravity in ADM formalism and constrained systems is described in more 
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detail in Refs. [45J. We want to point out that in this section the calculations are formal 
just as is the case of path integrals in field theory and in order to obtain some physical 
results additional considerations need to be implemented depending on the specific system. 



A. The Stratonovich-Weyl Quantizer 



Let F[hij, 7T*- 7 '; $, 71$ ] be a functional on the phase space T* (Wheeler's phase superspace) 
and let F[fi l \ Ay-; fi, A] be its Fourier transform. By analogy to the quantum mechanics case, 
we define the Weyl quantization rule as follows 



F = WiFihij,^ := fv V f^j V (J-^j V (£) F^ , Ay; /z, A]%* Ay-; A], 

(12) 

where the functional measures are given by Vhy = Yl x dhij(x) etc, {W|/t y ', Ay-; fM, A] : 
(/i JJ , Ay-; /i, A) G T*} is the family of unitary operators given by 

U\n % \ Xif /i, A] := exp |i J dx ^/^(x)hij(x) + Xij(x)^(x)+fx(x)^(x) + X(x)n^(x^^, (13) 

with hij, 7r u , $ and 7?$ being the field operators defined as 

hij(x)\hij, $) = hij(x)\h ij} $) and 7ry(ar)|7ry, 71$) = 7Ty(a;)|7ry, 71$), 

$(x)\hij,$) =$(x)\hij,$) and 7r*(a;)|7ry,7r$) = 7T*(a;)|7ry-,7r$) . (14) 

These states form basis satisfying the completeness relations while operators satisfy the 
usual commutation rules. Using the well known Campbell-Baker-Hausdorff formula, the 
completeness relations and the standard commutation rules we can write U\pL 1 ^ Ay-; /i, A] in 
the following form 



Xif fx, X] 



VhijVQ exp i^i J dxyjj 1 ' i: '(x)hij(x) + fi(x)^(x] 



hx 



hx 



ha + 



hx 



'j 



hx 



2 2 / \ J 2 2 

It is easy to show, from (15), that one can obtain the following properties: 



TrW[nV,X i fLi,X]} = 5[X ij }5 



Trj^V, X t f /i, x\-UW ij , K,; P*, A'] |> = *[A«,-Ay * 



2tt 

" ft 

2^ 



5[A]5 



2tt 



(15) 



(16) 



(/<'•' -/<' ; ') 



■5[A-A']5 



ft 

2^ 



(17) 



where Tr is the trace in some representation. 



Equations (15) and (12) lead to write 



(18) 



where the operator Q is given by 



x exp | — i J dxf/j, iJ '(x)hij(x) + Xij(x)n v '(x) + + \(x)ir<s,(x)^J X U\p? j , Ay-; A]. (19) 

It is evident that the operator is the Stratonovich-Weyl quantizer. One can easily check 
the following properties of 

(fi[^j,7T y ;$,7r«]) = (^[/i^,7r iJ ';$,7r$]), 



Tr{fi[^,7r^;$,7r $ ]} 



2vr/i 



- h'J ■ 5 



- 7T $ 

2tt/i 



<*[$-$']. 

(20) 



Now it is possible to obtain the inverse map of W by multiplying (12) by f2[/iy, 7r y ; $,7r$ 



taking the trace of both sides and using (17) we get 



F[hij, Tr"; $,tt*] = TW n[h ij} tt«; $, vr $ ]F k 



(21) 



One can also express fi[/iy, vr*- 7 ; $, 7r$] in a very useful form by inserting (15) into (19). 
Thus one gets 



ft[/^,7r^$,7r $ ] = / Vfa / 2^exp<j -- I dx[^(x)^(x)+^x)n iS> (x) 



X 



J 2 2/\ 3 2 2 



(22) 



B. The Star-Product 



Now we define the Moyal *— product. Let F = F[hij, ir^; $, 7r$] and G = G[hij, n lJ ; $, 7r$] 
be some functionals on T* that correspond to the field operators F and G respectively, i.e. 
F[hij,7r ij ; <I>,7r$] = W" 1 ^) = Tr(f2[/^-, tt«; $, tt $ ]F) and G^', 71$, $] = W^G) = 
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Tr(fl[7T w ', 71$, The functional corresponding to the product FG will be denoted by 

(F * G)[hij, ir^; $, 7r$] and after some long but direct calculations we have that 

(F*G)[/^,7r lJ ';$,7r $ ] := = Trjf2[/^, n ij ; $, tt*]FgJ, (23) 

gives rise to 

(F*G)[^,7r^';$,7r*] = F[^,7r^;$,7r*]exp|y ? 7r iJ ; $, vr,,], (24) 

where the operator P is defined as follows 



8hij(x) Sir^(x) 5iry(x)5hij(x)J J \8$(x) Stt^(x) 5tt^(x) 5&(x 

(25) 

This is precisely the Moyal *— product. 



C. The Wigner Functional 



In the deformation quantization formalism the Wigner function plays a very important 
role and in the same way finally we are going to define the Wigner functional. Let p be the 
density operator of a quantum state. The functional pw[h>ij, 7T y ; $,7r$] corresponding to p 
reads 

p w [hij, rr v ; $, 71$] = Tr jfif/i^, 7r %J ; $, 7r$]p 



6 



27r/l 



2vrn 



cxp 



C?X ( £ij(x)TT lj (x) + C,(x)lT^(x) 



(26) 



For a pure state of the system p — the equation ( |26| ) gives 



27r/l 



2nh 



exp 



I £ij(x)TT l:1 (x) + £(x)7T$(x] 



X ^* 



hi 



h-. + ^i- <i> . £ 

JJ + 2 ' 2 



(27) 



where $|^} = $] is the wave function of the universe. 

The expected value of an arbitrary operator F can then be obtained by means of p as 
Tr(pF) 



(F) 



Tr(p) 



(28) 
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It is now possible to write the equivalent of the constraints equations ([9]) in terms of the 
★—product and the Wigner functional as 



n±*p w [h ij ,'n ij iQ,Tr*\=0, (29) 

H i *p w [h ij ,ir i >]$,ir 9 ] = 7 (30) 
where H± and Hi are given by (J5f and Q respectively. We term Eq. (29) the Wheeler- 



DeWitt-Moyal equation. It is important to mention that the dynamics is completely gov- 
erned by these equations and that their explicit form depends of the particular system 
considered. These equations will be useful in the cosmological models which we will dis- 
cuss in the following section. Also, we want to note that the real and imaginary parts of 



Eq. (29), encoded in the *— product, are the generalized mass condition and the Wheeler- 



DeWitt-Vlasov transport equation presented in [2S] for flat minisuper space. 



IV. DEFORMATION QUANTIZATION IN THE MINISUPERSPACE 

Now we proceed to study the application of the deformation quantization construction 
to some models in quantum cosmology. The general construction exposed in the previous 
section in terms of functional integrals is defined in the whole superspace of 3-metrics and 
is necessary in a general case , however as a first approach we will deal with some models 
in the minisuperspace due to their simplicity and because they have been widely studied 
in the literature by other methods. Our intention is to obtain the Wigner function for 
these important cases and to motivate a further study using deformation quantization. We 
consider that this first step is necessary in order to gain some experience to eventually 
apply the deformation quantization formalism to curved phase-space. It is known that 
this formalism has a natural extension to these situations (see [221 123]) and allows another 
suitable extensions or generalizations. 



A. de Sitter cosmological model 

As a first example we will deal with a minisuperspace model where the phase space is 
bidimensional. We are going to calculate the Wigner function from the Wheeler-DeWitt- 
Moyal equation and also from its integral expression. Let start with a minisuperspace (where 
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the degrees of freedom are reduced to just one represented by the scale factor of the universe) 
corresponding to the Friedmann-Robertson- Walker metric 



ds 2 = l 2 p [-N 2 dt 2 + a 2 {t)dVtl] , (31) 

where a is the scale factor of the universe, N is the lapse function, dQ 2 is the metric of the 
unit three-sphere, l p = 2/3L p and L p denotes the Planck length. Introducing new variables 
q = a 2 and N = qN |54| the Hamiltonian H± cast out in an easy form 

n ± = \ (-4p 2 + Xq - 1) , (32) 

where A is the cosmological constant in Planck units. In the coordinate representation the 
Hamiltonian acquires a simple form and its dynamics corresponds to a particle in a linear 



potential. The Wheeler-DeWitt equation comes from the Hamiltonian constraint (10) and 
is written as 

4ft 2 ^ + Ag-l^(g)=0. (33) 
Depending of the boundary conditions it can be found [UJ the Vilenkin wave function 



* l ' W =2 Ai ((^J + 2 Bi ((2A«)575j- (34) 



the Hartle-Hawking wave function 



*™<»> =Ai {m^ ) • (35) 



and the Linde wave function 



*^ = - !Bl (»)' (36) 

where Ai(x) and Bi(x) are the Airy functions of first and second class respectively. 

One of the main differences between the Hartle-Hawking, Linde and Vilenkin wave func- 
tions is their behavior in the region q > 1/A. The Vilenkin tunneling wave function has the 
following expression 

*V = V-(<Z), (37) 

the Hartle-Hawking wave function is 

* HH = i/> + {q)+il>-(q), (38) 
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and the Linde wave function is written as 



(39) 



where ip-(q), i[)+(q) describe an expanding and contracting universe, respectively. In the 
classical allowed region q > 1/A the WKB solutions are 



tl>± = \p(q)] 1/2 ex P 



±i / p(q')dq' =F 



VK 



1/A 



(40) 



where p(a) = [— £/(a)] 1//2 = | [Ag — l] 1 / 2 . From the WKB wave function (40) we have 



pi/j±(q) « ±pip±{q). 



Taking into account that 



p = - 



4N' 



(41) 



(42) 



the Eq. (41) confirms the already given interpretation for ip±(q), i.e. negative values of 



p correspond to an expanding universe. Therefore, the Vilenkin wave function includes 
only an expanding component while the Hartle-Hawking and Linde wave functions include 
expanding and contracting universes with equal weight (for a different interpretation see 
[55]). In fact, for large values of q these wave functions have the following expressions 

(2A/1) 1 / 6 



2[vr 2 (Ag - l)} 1 / 4 



^HH{q) 



(2XH) 



1/6 



[^(Ag-l)] 1 /* 



cos S 



(2Xh) 



-iS 



1/6 



2[7r 2 (Ag- 1)]V4 



(e lS + e-* s ) 



*L{q)ni r S UH)1/ L .. sin S 



where 



[vr 2 (Ag- 1)]V4 
5 



2[vr 2 (Ag- 1)] ] 



(43) 

(44) 
(45) 

(46) 



Now we proceed to calculate the Wigner function directly by solving the Wheeler-DeWitt- 
Moyal equation (29) 

>H ± *pw = 0- (47) 
where the corresponding Moyal *— product is given by 



'iH « 
■k = exp { — V 



exp 



ih I d d 3d 
2 1 dq dp dp dq 



(48) 
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Taking in consideration the exponential power series expansion we have 

1 



-Ap z + \q-l) 



OO j 



n=0 



/?■! 



PW 







which can be written as 



ih 



Ap + \q- l)p w + —XdpPw + AiHpd q p w + h d p w = 0. 



(49) 



(50) 



If we define a new variable z = Ap 2 — Xq + 1 the imaginary part of the former equation is 



identically zero and for the real part we obtain a new form of the equation (50) 

2 X 2^PW n 
Z/0H/ = ' 

whose solution is 



p w = ciAi 



Ap 2 - Xq + 1 



(51) 



(52) 



In order to obtain the last result the formalism assumes the existence of the Wigner trans- 
form, where the range of integration in this transform is from minus to plus infinity in 
the variable q. However, the valid range interval for the variable q is positive but still the 



Wigner function is very similar to (52) because of the exponential decay for positive val- 



ues of the argument of Ai(x). Last equation indeed admits another solution corresponding 
to Bi ^ ^^'2/3" 1 J ? Du t in this case, the Airy function Bi(x) is divergent for positive x and 
this part cannot be included as a suitable Wigner function. In fact, in order to obtain an 



appropriate solution from (47), the potential could be complemented by an infinite barrier 



avoiding the existence of negatives values for q. However the implementation of this proce- 
dure is very cumbersome as it has been shown in [56] and [57] and it is not convenient to 
develop it here. 

Instead of this we can calculate the Wigner function using the following integral expression 



Pw(q,p) 



27iH 



exp 



(53) 



Employing the convolution theorem and the Fourier transform of the Airy function we get 
the Wigner function for the Hartle-Hawking wave function 



Pw{q,p) 



2 l/3 



rAi 



Ap 2 - Xq + 1 

(frA) 2 /3 



(54) 



7i(nxy/ 3 ' 

The last result was obtained integrating out from minus to plus infinity in the variable q (in 
fact this result was already derived in [30J ) . However, the wave function is only valid defined 
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FIG. 1: The Hartle-Hawking Wi gner function (ft = A = 1). FIG. 2: Hartle-Hawking Wigner function density projec- 
The figure shows many oscillations due to the interference be- tion. It can be observed that the classical trajectory does not 
tween wave functions of expanding and contracting universes. coincide with the highest peaks of the Wigner function. 

for positive values of the scale factor and the q variable. Due to the fact that the Airy 
function Ai(x) presents an exponential decay for positive values of x, the Wigner function 
is very similar to the expression obtained in terms of the Airy function. If we restrict the 
range of integration of q for positive values the computation for an analytical expression of 
the Wigner function is very complicated. This problem is similar as the one we faced for 
the Wheeler-DeWitt-Moyal equation. We cope with this complication by implementing a 
Fortran code in order to calculate numerically the Wigner function. The result is given in 
Fig. [I] and Fig. [2j where the continuous line corresponds to the classical trajectory (56) and 



the dashed line to the trajectory (57) 



For the Vilenkin wave function it is very difficult to get an analytical expression for 



the Wigner function directly from (53) even using the WKB approximation for the wave 
function taking into account that the integration should be restricted for positive values of 
q. Again, we performed a numerical analysis to obtain the results that are depicted in Fig. 
[3] and Fig. [4j where the meaning of the continuous and dashed lines are the same as in the 
Hartle-Hawking case. 

As before, the Linde Wigner function is very complicated to work with. We calculated 
the Wigner function numerically and its behavior is presented in Fig. [5] and Fig. |6j 

These plots can be described using the WKB approximation for the wave function and 
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FIG. 3: The Vilenkin Wigner function. It is observed a FIG. 4: The density projection of the Vilenkin Wigner 
clear maximum and less oscillations compared with the Hartle- function. The classical trajectory is at some parts on the 
Hawking case. maxima of the Wigner function and has only one branch. 




P 




FIG. 5: Wigner function for the Linde wave function. The FIG. 6: The density projection of the Linde Wigner func- 
figure shows a reduction in the amplitude of the oscillations tion. In this case the classical trajectory coincides with the 
compared to the Hartle-Hawking. highest peak of its corresponding Wigner function. 

the general statements presented in [30]. In terms of the WKB density matrix Pwkb the 
evaluation of 

1 f 00 

Pw(q,p) ~ — r / daexp(-2ipa)p WK B(q + a,q - a) (55) 

™ J-oo 

for the Hartle-Hawking wave function includes four terms of the form: exp(iS(q + a) — iS(q — 
a)), exp(—iS(q+a)—iS(q—a)), exp(iS(q+a)+iS(q—a)), exp(—iS(q+a)+iS(q—a)). These 
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terms come from the products ip + (q + a)ip+(q — a), ip-{q + a)ip* + {q — a), ip + {q + a)ip*_ (q — a), 
ip-(q + a)i}j*_(q — a), respectively. By means of the stationary phase approximation, each 
one of the former terms will contribute to the Wigner function in different regions bounded 
by the classical trajectory 

P 2 = \{M-1\ (56) 

and by the trajectory 

P 2 = \{M-1\ (57) 

which corresponds to points where the path integral has a saddle point at zero momentum 
[50] . i.e. where the classical action changes sign. In the upper region of the phase space 
between the classical trajectory and the dashed curve the Wigner function has contributions 
mainly coming from the first term in the density matrix. The region inside the dashed curve 
receives the principal contributions form the second and third terms of the density matrix. 
The fourth term of the density matrix has a saddle point in the bottom part of phase space 
between the classical trajectory and the dashed curve and as a consequence it contributes 
predominantly in this region (see Fig. [2J. 

The Linde Wigner function has a similar structure like the Hartle-Hawking case but 
there is one important difference: the interference terms between contracting and expanding 
universes have an opposite sign with respect to Hartle-Hawking. The consequence of this 
difference is the reduction of the amplitudes of the oscillations inside the classical trajectory 
(see Fig. [5J. 

The Wigner function for the Vilenkin wave function only receives a contribution from 
density matrix corresponding to tp-(q + a)tp*L (q — a) and it has a predominantly contribution 
only in the lower region of phase space between the classical trajectory and the dashed line. 
This behavior is reflected in Fig. |4| It is important to mention that decoherence of the 
Vilenkin Wigner function appears to be easier taking into account that large amplitude 
terms due to interference between collapsing and expanding universes are not present near 
p = 0. 

We discuss now the interpretation of these numerical results. For the Hartle-Hawking 
Wigner function we see an oscillatory behavior where the largest peaks are not on the 
classical trajectory but away from it by a distance of O^ 2 / 3 - 1 ), we can appreciate this fact 
in the density plot. It should be remarked the existence of oscillations of the Wigner function 
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that are not near of the classical trajectory, however the heights of these peaks decreases 
with the distance to the classical trajectory. 

The Wigner function for the Linde wave function present a similar pattern, in general 
terms, like the Hartle-Hawking wave function, however it is possible to see some differences. 
We can appreciate more fluctuations of the Wigner function inside the region corresponding 
to the classical trajectory than the Hartle-Hawking Wigner function. Furthermore, the 
amplitude of the oscillations are smaller for the Linde wave function than for Hartle-Hawking. 
The largest peaks of the Wigner function correspond to the classical trajectory. 

In the case of the Wigner function for the Vilenkin wave function the classical trajectory 
is at some parts of the curve at the middle of the largest peaks of the Wigner function, i.e. 
the Vilenkin Wigner function has the largest peaks more closely to the classical trajectory 
than the Hartle-Hawking Wigner function, but only in one part. This behavior is expected, 
as we explained above, due to the fact that the Vilenkin wave function represents the tun- 
neling wave function and as a consequence there is only one part of the classical trajectory 
corresponding to negative values of the momenta (expanding universe). 

B. Kantowski-Sachs model 

Another interesting case to deal with under the deformation quantization procedure is 
the cosmological Kantowski-Sachs model (KS) [HE]- We consider the metric in the Misner 
parametrization [59] 

ds i = _jv 2 dt 2 + e 2 ^dr 2 + e~ 2V ~ 3l3 e- 2V ~ 3n (d9 2 + sin 2 9d<p 2 ). (58) 

Choosing a particular factor ordering the Wheeler-DeWitt equation takes the following form 

(-P^ + Pj - 48 exp(-2v / 3fi)) *(0, /?) = 0, (59) 

where Pq = —ih-^ and Pp = —ih^. The solutions of the former equation are given by [59] 

*±(0, P) = e ±l ^K ifi Qe-^) , (60) 

where Ki V [x) is the MacDonald function of imaginary order. The normalized wave function 
is 

W) ^f-^) Kii ^) (61) 
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which satisfy [60 



dndf3^i(n, (3)y u ,(n, (3) = 5{u 2 - v' 



(62) 



In order to write the Wheeler-DeWitt-Moyal equation (29) it is very useful to employ the 
next relation 

f(x,p)*g(x,p) = f ( x+ —d p ,p- -^d x j g(x,p). (63) 
In this way, we obtain the following equation 



P n -p,)\(p,-pX -^<" + ^ V/ MaFn ,^ e) = , (64) 



which can be split into two equations corresponding to its real part 

h* „ „. h 2 



-P 2 + \dl + P 2 - \dl - 48e- 2 ^ cos (V3W Pa 
and its imaginary part 



p = 0, 



(65) 



KPndn) ~ KPpdp) + A8e~ 2V ~ 3Q sin (V3hd P ^ 



p = 0. 



(66) 



If we propose p(Q,Pn,(3,Pp 



Pn(Q, Pn)p/3((3, Pp), and taking into account that 



f(x) = f(x + i\/3h) and also that for free particle in the f3 parameter dppp = 0, 



we get from equation ( 66 ) 



48\/3? 



hPn 



e - 2 ^( P (fi, p n + iV3h) - p(n, p n - iVzti)) 



576 e- 2 ^ n T 



P n + 2iVSh) - 2p(Q, P n ) + P n - 2iV3H) 
^— (p(Q, P Q + 2iV3h) - p{n, Pn - 2iV3H)) 



(67) 



Using the last equation in (|65|) we finally obtain 

pn(n, P n + iVSH) - pn(n, P n - iVSH) 



o l2VSih 2 e- 2 ^ n 



+ 



PnPn 

[Pj + 3^ 2 ) 
Pn(P^ + 3h 2 ) 



Pn 

p u (n, P n + 2iV3h) - 2p n {ft, Pn) + Pn(n, Pn - 2i\f?>h) 
pn(n, Pn + 2iVSH) - (^(£1, P n - 2iy/3h) 



- 24e 



pn(fi, Pn + iV3h) + pa(n, P n - iV3h) 



-P 2 



pPn- 



(68) 
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It is hard to obtain directly a solution to this equation, so we will follow a different approach 
and will use the integral representation for the Wigner function. Using the following result 
(see Sec. 19.6 formula (25) in [51] and the comment in [5"2"] ) 



dwiwzY^w^K^a/w)^ 



wz 



2-°-^ 2 a a Gf A 



a 2 z 2 
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[L — a —fi — a I v 1 v 
2 ' 2 ' 4 + 2' 4 ~ 2 



where G$ ( ^ 



(69) 



/J — IT — /-t— U III' 1 f 



. ,, . 4 + ^, | — | ] is a special case of Meijer's G function (see Sec. 

5.3 in [63]) 

a i; z = 1, 

6j, j = !> -,3 
we calculate the Wigner function for the Q part 



(70) 



3 1 / 2 , /7TI/ 



ft 



(71) 



Then, we obtain the following expression for the Wigner function 



7T" 



16fr 2 l^ 2 



10 / 16 _4^30 



X °ol I 



1 (-v Pn\ 1 zi/ 1 



4 ' V 2 ^ \/3^7 '4 ' " V 2 ^ ' V3^7 '4 ' 2ft' 4 2^; 



(72) 



Now, employing the Meijer's function property 



o», i = 1, -,P 



bj, j = 1, 



cii + a, i = 1, 
+ CT > J = l,-,5 



(73) 



the equation ( 72 ) can be written in the following form 



Pn{n,P n ) = 



40 I 16„-4\/3n 



x G ol \ ^ 



+ i 



Wh 2 
v 



Pn \ 1 , . / -/' 



2h 2y/3hJ ' 4 



Pn \ 1 
2fi + 2y/3h ) ' 4 



Pn \ 1 



(74) 



It is possible to verify that the Wigner function indeed satisfy equation (68) 



In order to extract physical information we plot the Wigner function for several values 
of v. We can say from Fig. [9] and Fig. 10 that the classical trajectory is near the highest 
peaks of the Wigner function for values close to v = 1. For values of v smaller than one we 
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FIG. 7: The Kantowski-Sachs Wi gner function for u = 0.5. FIG. 8: Kantowski-Sachs Wigner function density projec- 
Few oscillations are present for this case with a clear maximum tion for u = 0.5. It can be observed that the value of the 
around f2 = 6 and Pq = 0. Wigner function is large in an ample area. 

can see from Fig. [7] and Fig. [8] that there are less oscillations but the classical trajectory 
does not correspond to the highest peaks, in fact, there is an ample region where the Wigner 



function is large. For values of v bigger than one the Fig. [TI] and Fig. [12] show an increment 
in the number of oscillations of Wigner function and the peaks of the oscillations are far 
away from the classical trajectory. We can conclude that the quantum interference effects 
are enhanced for larger values of v. 

Recently the noncommutative version of the Kantowski-Sachs model have been studied in 
[21] and it turns out interesting to analyze it using the deformation quantization formalism. 
Now from the previous results we can present the noncommutative version of the Kantowski- 
Sachs model. We consider a Kantowski-Sachs cosmology with a noncommutative parameter 
9. The operator algebra in the phase space is given by 

[z a , n Za ] = iHS aj3 , 

[Z a ,Z p ]=i9e a p, [n^,ng =0, (75) 

where Z a = (f2, (3) and I1^ Q = (Pf^P^). Of course further generalizations can be imple- 
mented in the general case with [II^ a , Hz*] 7^ 0. 

The noncommutative Wheeler-DeWitt equation is written as 

' +!^ + 48exp(-2v / 3ft + v / 3#P (3 )) = 0, (76) 



on 2 d(3 2 
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FIG. 9: The Wigner function for the Kantowski-Sachs wave FIG. 10: Density projection of the Kantowski-Sachs 
function for u = 1. The number of oscillations are maxima Wigner function for u = 1. It can be observed that classi- 
around Pq = 0. cal trajectory is close to the exterior peaks of the oscillations. 

3 




Pn 




1 2 3 4 5 6 7 



FIG. 11: Kantowski-Sachs Wigner function for v = 4. FIG. 12: Kantowski-Sachs Wigner function density pro- 
There is a considerable increase in the number of oscillations jection for v = 4. The classical trajectory (near the Pq, axis 
centered at Pq = 0. at the upper and bottom part of the Fig.) is far away from 

the peaks of the Wigner function. 



where P, 
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-ifr-§p- The solutions are given by 



^(n,/3) = e ± ^K q l -exp 



(77) 
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The normalized wave function is 



(78) 



which satisfy 



dftd/3tf * (ft, /^(ft, /?) = 5(u 2 - v' 1 ). 



(79) 



Following a similar procedure as for the commutative case we are able to write the Wheeler- 
DeWitt-Moyal equation in the following form 



— 48 exp 



(80) 



and the corresponding difference equation is 
9 12V3th 2 e- 2 ^( n ^ e ) 



Pn 



P 9 n (n, P n + iV3h) - p e n (n, P n - iVsh) 



+ 



(pi + m 

Pn(P£ + 3h 2 ) 



+ 2i>/3n) - 2p^(fi, P n ) + Pn - 2i<s/3fi) 



Pn + 2iy/Zh) - P e n (n, P n - 2iV?>h) 



= -Pi 



2„0 



(81) 



The solution to this equation is again difficult to obtain in a direct way but employing the 
integral representation of the Wigner function we can calculate the corresponding noncom- 
mutative Wigner function from 

(82) 



Thus we find 



7T> 



16^ 2 



Tl 2 



2iP n 

V3(f}=F^0)\ ^ 



x GS2 ( ^e-^K^^) 



1,./" , p ^ "\ 1 



4 + ^ + ^'4 + *'-- H - 1 



^ Pn \ 1 ii/ 1 if 



(83) 
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FIG. 13: The Kantowski-Sachs Wigner function weighted FIG. 14: Kantowski-Sachs Wigner function density pro- 
by a Gaussian function. It can be observed the Gaussian form jection weighted by a Gaussian. There is a clear maximum 
between Q = 1 and Q = 2 centered at Pq = with some very with some other transversal oscillations around Pq = 0. 
smooth oscillations around Pn = 0. 

Just as in the previous commutative KS case we can rewrite the former equation as 

ihyKV/h) e V T 2 ) 



Pn {U, Pa) 



Wh 2 

X fv Pg \ X 
4 \2h 2V3h) ' 4 



i , Pn \ X . / v Pu \ 1 



— v 

2h ' 2-ShJ ' 4 



-if _ iPn \ 
2 v / 3/lJ 



(84) 



We can verify that this expression fulfils the equation (81) and the noncommutative effect 
corresponds to a displacement in the Q variable. 

The noncommutativity consequences can be seen more transparently if we consider a 
wave packet weighted by a Gaussian 



r 



(85) 



The resulting values of Eq. (85) are introduced in the Wigner function (53) and the 



integral is performed numerically for a value of (3 = 0. The result is plotted in Figs. 15 and 



16 for the values a = 4 and 6 = 1. Figs. 13 and 14 correspond to 9 = and Figs. 15 and 
16 correspond to 9 = 4. Here we are interested in describing which is the effect of the 9 



parameter. Figs. 13 and 14 shows that for 9 = there is only one preferred semi-classical 



state of the universe. Figs. [15] and [16] show a significative difference consisting in having 
several smaller picks representing semiclassical states of lower probability. The additional 
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FIG. 15: Kantowski-Sachs noncommutative Wigner func- FIG. 16: The density projection of the Kantowski-Sachs 
tion weighted by a Gaussian. Several peaks of different am- noncommutative Wigner function weighted by a Gaussian, 
plitudes in the Wigner function can be observed around the Due to the noncommutativity more peaks and oscillations ap- 
Gaussian form. These are interpreted as different universes pears around Pq = 0. 
connected by tunneling. 

states contribute to a landscape of vacua since these new vacua can be reached by tunneling. 
Summarizing, at the semi-classical level, it is confirmed the observation made in Ref. [21] 
concerning that the noncommutativity of the minisuperspace leads to new possible states of 
the universe which contribute to its evolution at early stages. 

C. String cosmology with dilaton exponential potential 

Now we can treat the case of string cosmology (HI] with dilaton exponential potential in 
the same way as in the previous subsection. We consider the D = 4 tree level string effective 
action 

S = J d'x^-g [n g - l -(d<P) 2 - V{4>)\ , (86) 

where is the Einstein frame metric, g is its determinant, TZ g is the Ricci scalar compatible 
with g^, (ft denotes the dilaton and we consider the dilaton potential V(<p) = V^e^. We use 
the following metric 

ds 2 = ~^^dt 2 + a[t)% j dx i dx i 1 (87) 
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and N(t) is the lapse function. By means of the following new variables 



x = - In 
6 



(u + v ){u — v)' 



z = — In 
6 



2(u + v) 
(u — v) 2 



(88) 



where u = f^cosh20j and v = f^sinh20Y the Wheeler-DeWitt equation is written 
the gauge N' 1 = ^(w — t>) as follows jHS] 

d 2 d 2 



m 



H±ip(x, z) 



0. 



(89) 



The solution is expressed in terms of Bessel function in the x variable and as a free wave in 
the z variable, therefore the wave function is 



iP(x,z) = e ±lkz J ±lk/3 WV e* x ) 



(90) 



where k is the separation constant. We choose the delta function normalization for the wave 
function. In this case it is convenient to use the wave function in the x variable 



(91) 



2y sinh(7r/c) 

where it satisfies (see section (4.14) in [66] and [6j 

dxr(x) ±k *lj(x) ±k , = 5(k 2 - k' 2 ). (92) 
We proceed now to obtain the expression for the Wheeler-DeWitt-Moyal equation, again 



it is convenient to use (63). We get 
ih 



P: 



d z ) + 9V e 



which can be separated in its real part 



-Pi + + Pl - ^d 2 z + 9V e«* cos (3/tfp.) 



p(x,P x ,z,P z ) = 0, 



p(x, P x , z, P z ) = 0, 



and its imaginary part 

[h(P x d x ) - h(P z d z ) + 9V e 6x sin (3hd Px )] p(x, P x , z, P z ) = 0. 



(93) 



(94) 



(95) 



If we propose p = p x (x, P x )p z (z, P z ), and consider that for a free particle in the z parameter 



d z p z = 0, we get from equation (95) 

27 

hP x d 2 xPx (x, P x ) = - —iV e &x { Px {x, P x + 3ih) - p x (x, P x - 3ih)) 



81 {V e 



6x\ 2 



P x (p x (x, P x + 6ih) - 2p x (x, P x ) + p x (x, P x - 6ih)) 



Ah (P 2 + 9h 2 ) 

+ 3iH(p x (x, P x - 6ih) - p x (x, P x + 6ih)) 



(96) 
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Using the former equation in (94) we finally obtain 



-P*p x (x, P x ) = - P 2 x p x {x, P x ) + 



27ihV 
4R 



e 6x (p x (x, P x + 3ih) - p x (x, P x - 3ih)) 



+ ^^(p^x, P x + 3th) + Px (x, P x - 3th)) 
(9V e 6x )' 



zt^\(Px — 3ih)p x (x, P x + 6ih) 
4 2 P X (P 2 + 9h 2 ) lK )Hx\ , x ) 

2P xPx {x, P x ) + (P x + 3th)p x {x, P x - 6th)}. 



(97) 



This difference equation is complicated to solve in a direct way so as with the Kantowski- 
Sachs case we use the integral representation of the Wigner function. Employing the follow- 
ing result (see Sec. 19.3 formula (45) in [61]) we can calculate the Wigner function in terms 
of the Meijer's function 



x^J^^Jvibx-^dx = 2 p - 1 a- p G™ 
Then we obtain, 
Px(x,P x ) 



a 2 b 2 



16 



v p + p p — p v 
2' 2 ' 2 '~2 



(98) 



8hiT sin(7r/c) 

/T/2 12a: 

± G, 20 ' 



T/2 12a; 
n 20 I v e 



r 01 



16fr 4 



l -(-k + P x ), '-(k - P x ), % -{-k - P x ), % -{k + P X/ 



r 04 



16fr 4 



'T/2 12a; 

± G?° 1 



r 04 



r 



16fr 4 



T/2 12a; 
20 / v e 
04 



16/l 4 



l -{-k + P x ), % ~{-k - P x ), % -{k - P x ), % -(k + P x ] 
l -{k + P x ), '-(k - P x ), '-(-k - P x ), '-(-k + P X/ 
l -(k + P x ), j(-k - P x ), % -{k- P x ), % -{-k + P x ] 



(99) 



We want to note, that it is also possible to write the Wigner function in terms of the 
hypergeometric function F 3 employing [58] ■ This is an straightforward but long calculation 
that we will not present here. 



It can be verified that p x given by Eq. (99) satisfies the equation (97). 



D. Baby Universes 

In this subsection we will consider another example of the use of the deformation quan- 
tization formalism to wormhole solutions in general relativity [69TT7T] . Here we will have a 
system with two coordinates of the flat minisuperspace (2 degrees of freedom). 

We are going to consider a baby universe with conformal matter 0o and a three-metric 
hij defined on a Cauchy hypersurface S of a closed wormhole universe. The matter is 
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represented by a conformal invariant scalar field expanded in hyper-spherical harmonics Q n 
of the surface S given by 4> = - J2 n fnQni where a is the scale factor and f n are orthonormal 
modes. The metric on S is given by h{j = a 2 ■ (O^- + e^). Here fiy is the metric of a unit 
three-sphere S 3 and = J2 n ( a rfoijQn + KPijn + c n S ijn + d n G ijn ) . The Q„'s are the scalar 
harmonics of the 3-sphere, Pij n is a suitable combination of Q n , Sij n is defined in terms of 
the transverse vector harmonics and Gij n are the transverse traceless tensor harmonics on 
S. 

On the gravitational part, in a suitable gauge (a n = b n = c n = 0) of on S and 
considering the case without gravitons (d n = 0), one can express the three metric simply as: 
hij = a 2 ■ flij. Thus, the wave function \1/ is a function of the scale factor a and the harmonic 



modes of the scalar field f n . This wave function fulfills the Wheeler-DeWitt equation (10) 
of the form 



£(-S + " 2/ »)-(-S +a2 



*(o,/ n ) = 0, (100) 



where we have implemented the canonical relation on the momenta pf n = —ifi-^- and 

In the context of quantum cosmology the solution factorizes in a purely gravitational 
part and an purely matter part. Both of them correspond to harmonic oscillators and the 
solution is given by 

*(a,/ n ) = AH m (a)ex V f-^pj ■ J] H mn {fnM exp [~\ , (101) 

where H m (x) are the Hermite polynomials and A is a normalization constant. This solution 
satisfies the boundary conditions such that: ip(a, f n ) — » as a — > oo and it is regular at 
a = 0. 

Here it is assumed, as in |70j, that the zero-point energy of the gravitational sector 
will precisely compensates the zero-point energy of the matter oscillators as it happens in a 
supersymmetric theory. This solution represents a closed universe carrying m scalar particles 
in the n-th mode. Thus the ground state |^o) wn l correspond to m = and n = 0, i.e., the 
absence of matter particles and consequently excited states of the scale factor part. 

The wave function factorization \P(a, f n ) = ipoi^'Yln ^nifn) comes from the inner product 
between (a, f n \ and \ip , ip n ) . The ground state is given by |\J/ ) = \i J o)a ® IV'o}/- 
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Let b and the annihilation and creation operators of the m a -th modes of the gravi- 
tational sector defined by: 6|0 a ,m/) = and [6^] ma |\l/o) = |^a; m /}- Now, let d and d) 
the annihilation and creation operators of the m/-th modes of scalar particles defined by: 
d\m a ,0f) = and [<v] m f \^ ) = \m aj mf). The combination yields b- d\0 a ,0f) = and 
p} m - ■ [d}} m f\^ ) = \m a ,m f ). 

In the WWGM formalism we have the following general stationary Wheeler-DeWitt- 
Moyal equation 

H B *Pw(a,Pa,fn,Pf n ) = 0, (102) 

where pw = pw(a,Pa, fn,Pf n ) is the Wigner function and Hb stands for the baby universe 
Hamiltonian 

H B = H f + H a , (103) 

where 

H f = E (Pi + "nfn) ( 104 ) 
n 

with u n = n and 

H a = p 2 a + a 2 . (105) 

The Moyal product is given by 

f*g = fexp( l ^VB)g ) (106) 



v 2 

where the corresponding Poisson operator has the following form 

t> Jp J r %=X^(A I I <9 \ + d_Jl d_d_ 

3 f „\dfndp fn dpf n df n J da dp a dp a da 



Therefore we have that the equation (102) at the order h can be written as 



E (<«*/» 



(107) 



We can write down p wo = p^ ■ p{y then we can separate (102) into two parts 

Ha* Pwo(^Pa) = -Ep a wo (a,p a ), (108) 

and 

H, * p' wa (f,P,) = Ep{ vo (f,p f ). (109) 



9pL ^)_ a .^fa. +fc .$!&s = 0i (110) 



dPU n dfn J dp a da 
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where pwo is the Wigner function for the ground state. Thus the solution to this equation 
is given by 



PW0(Pa,a,Pf n ,fn) = Pwo( a 'Pa)- Pwo(fn,Pf n ) =Aexp 

The density matrix for all excited states is given by 



lpl + a 2 ) 



n 



exp 



Pr,s 



P] r .[?]'|*o><*o|[^]'-P t ] r . 



(Ill) 



(112) 



The WWGM formalism allows to compute from this density matrix the general Wigner 
function for all excited states 



[pw]m,n = [&*] * • • • * [&*] * [<f ] * • • • * [<T] * p wo * [<T] * • • • * [d*] * [b*] * • • • * \b*]. (113) 



It is straightforward to show that it leads to the solution 



Pw {Pa,a, f n )Pfn ) = A exp 



(114) 

Here L m (x) is the Laguerre polynomial of degree m. Remember that they are 
related to the Hermite polynomials through the familiar formula L n (x 2 + y 2 ) = 
(-l) n 2 _2n YZi=o m \{n- m )\ H2m{.x)H 2 n-2m{y)- The case for the minimal coupling scalar matter 
follows a similar treatment and will not be discussed here, but analogous formulas can be 
obtained for this case. 



V. FINAL REMARKS 

In this paper we have constructed the WWGM formalism in the flat superspace (and phase 
superspace). The WWGM correspondence is explicitly developed and the Stratonovich-Weyl 
quantizer, the star product and the Wigner functional are obtained. These results can be 
used in general situations but in a first approach we applied the formalism to some interesting 
minisuperspace models widely studied in the literature, in particular, we used the Moyal star 
product to describe some relevant cosmological models in phase space. 

We studied de Sitter quantum cosmology using the Hartle-Hawking, Vilenkin and Linde 
boundary conditions, where we have found numerically the Wigner function for all of these 
cases. 
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For the Hartle-Hawking wave function the behavior of the Wigner function presents 
many oscillations due to interference terms between the wave functions of expanding and 
contracting universes. Similarly as in Ref. [30] our result shows that the highest peaks of 
the Wigner function do not coincide with the classical trajectory of the universe. 

The Linde Wigner function has a similar behavior to the Hartle-Hawking case, where there 
are also expanding and contracting components of the wave function. The main difference 
in their corresponding Wigner functions is just a sign in the interference terms between the 
expanding and contracting wave functions and as a consequence it produces a reduction in 
the amplitude of the oscillations inside the classical region. 

In the case of the Vilenkin tunneling wave function we notice that the Wigner function 
has less oscillations compared to the Hartle-Hawking Wigner function. This is explained by 
the fact that there are less oscillation effects because there is only an outgoing wave. In this 
case the classical trajectory corresponds almost to the maxima of the peaks of the Wigner 
function. 

However the classical limit for these three models is difficult to obtain due to the existence 
of oscillations in the Wigner functions. It is important to note that decoherence of the 
Vilenkin Wigner function is in principle simpler to obtain since the interference terms are 
absent because there is only an expanding universe around p = 0. 

For the Kantowski-Sachs cosmological model we found the Wheeler-DeWitt-Moyal equa- 
tion which is equivalent to a differential-difference equation. We found its exact solution in 
terms of the Meijer's function. We observe that the classical trajectory corresponds to the 
highest peaks of the Wigner function for values close to v = 1. The situation for v < 1 in the 
Wigner function corresponds to have less oscillations and the classical trajectory does not 
correspond to the highest peaks. For values of v > 1 we have an increment in the number 
of oscillations of the Wigner function and we do not have the peaks of the oscillations near 
the classical trajectory. Thus it seems that v could be regarded as a parameter encoding 
the quantum interference effects. 

We have also considered the noncommutative Kantowski-Sachs model. In a similar way 
we obtain the analytic Wigner function and its differential-difference equation. This case 
presents a noncommutative parameter 9 deforming the Wheeler-DeWitt-Moyal equation. 
The Wigner function is determined in terms of the above Meijer's function with shifted ar- 
gument by a factor proportional to 9. We have constructed numerically the Wigner function 
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with wave packet weighted by a Gaussian to see the effects of the noncommutativity. There 
were found several peaks in the Wigner function around the Gaussian which can be inter- 
preted as different universes connected by tunneling. Thus, at the semi-classical level, the 
statement made in [21 J about that the noncommutativity of the minisuperspace produces 
new possible states of the universe at early stages is confirmed. 

String cosmology with dilaton exponential potential is also discussed. We found the cor- 
responding Wheeler-DeWitt-Moyal equation and the equivalent differential-difference equa- 
tion. These equations have an exact solution in terms of the hypergeometric and Meijer's 
functions. 

Baby universe solutions are also obtained in this context where the Wigner function is 
calculated by finding the exact solution of its Wheeler-DeWitt-Moyal equation consisting in 
two decoupled deformed harmonic oscillators in terms of Laguerre polynomials. 

It is important to remark that this work opens the possibility of treating several im- 
portant questions that remain unsolved in quantum cosmology with a novel approach. For 
instance, deformation quantization allows to deal with systems having phase spaces with 
nontrivial topology. In quantum cosmology, the existence of symmetries implies that the 
phase-superspace and in particular, the phase-minisuperspace will be reduced by the im- 
plementation of these symmetries leading to a non-trivial topological space with a non-flat 
metric. Therefore, deformation quantization is able to treat these mentioned cases in a nat- 
ural way. Moreover the extension to supersymmetric quantum cosmology [TTl EH [20] can 
be also treated applying the results of 



Another point to remark is that the Wheeler-DeWitt-Moyal equation (29) proposed in 
the present paper contains the generalized mass-shell equation and the Wheeler-DeWitt- 
Vlasov transport equation in the flat minisuperspace [29] encoded in its real an imaginary 
parts respectively. The former result was obtained using the Moyal *— product and then all 
the technics of deformation quantization developed for a long time can be apply to it. It 
is known that the mass-shell equation and the Wheeler-DeWitt-Vlasov transport equation 
admits a suitable generalization to curved minisuperspace [29]. It would be interesting to 
deal with curved symplectic cases where the Fedosov's approach could be applied in order 
to find the Wheeler-DeWitt-Moyal equation in non-trivial phase superspaces and obtain 
the mass-shell and the Wheeler-DeWitt-Vlasov equations. We will study this problem in a 
further communication. 
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Besides, as was mentioned before, the problem of obtaining the classical limit by imple- 
menting a coarse graining is relevant. It is possible to model a coarse graining through the 
Liouville equation with friction and diffusion terms [30], this approach can be addressed in 
the context of the deformation quantization formalism. 

To conclude, we consider that deformation quantization possesses various advantages 
in order to deal with more complicated problems in quantum cosmology for example, to 
treat systems with non trivial topology or with curved minisuperspaces. For these cases 
the canonical quantization could lead to the existence of nonhermitian operators which is 
avoided in deformation quantization as a result of the use of classical objects instead of 
operators. For these reasons more examples and further research is needed to develop this 
approximation. 

Acknowledgments 

The work of R. C, H. G.-C. and F. J. T. was partially supported by SNI-Mexico, CONA- 
CyT research grants: J1-60621-I, 103478 and 128761. In addition R. C. and F. J. T. were 
partially supported by COFAA-IPN and by SIP-IPN grants 20100684 and 20100887. We 
are indebted to Hector Uriarte for all his help in the elaboration of the figures presented in 
the paper. 



[1] B. S. DeWitt, Phys. Rev. 160, 1113 (1967). 

[2] E.P. Tryon, Nature (London) 246, 396 (1973). 

[3] P. I. Fomin, Dokl. Akad. Nauk Ukr. SSR 9A, 831 (1975). 

[4] D. Atkatz and H. Pagels, Phys. Rev. D 25, 2065 (1982). 

[5] A. Vilenkin, Phys. Lett. B 117, 25 (1982). 

[6] L. P. Grishchuk and Ya. B. Zel'dovich, in Quantum Structure of Space and Time edited by 

M. Duff and C. Isham, (Cambridge University Press, Cambridge, 1982). 

[7] J. B. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 (1983). 

[8] A. D. Linde, Lett. Nuovo Cimento 39, 401 (1984). 

[9] V. A. Rubakov, Phys. Lett. B 148, 280 (1984). 

33 



[10] A. Vilenkin, Phys. Rev. D 30, 509 (1984). 
[11] A. Vilenkin, Phys. Rev. D 50, 2581 (1994). 

[12] R. Bousso and S. W. Hawking, Phys. Rev. D 54, 6312 (1996); J. Garriga and A. Vilenkin, 
Phys. Rev. D 56, 2464 (1997); A. D. Linde, Phys. Rev. D 58, 083514, (1998); N. G. Turok, 
S. W. Hawking, Phys. Lett. 432, 271 (1998); A. Vilenkin, Phys. Rev. D 58, 067301 (1998). 

[13] R. Brustein and S. P. de Alwis, Phys. Rev. D 73, 046009 (2006). |arXiv:hep-th/0511093| . 

[14] C. Kounnas, N. Toumbas and J. Troost, JHEP 0708, 018 (2007). |arXiv:0704.l"996| [hep-th]]. 

[15] M. Gasperini, Elements of string cosmology, Cambridge, UK: Cambridge University. Press. 
(2007). 

[16] A. Strominger, in Quantum Cosmology and Baby Universes, eds. S. Coleman, J. B. Hartle, 
T. Piran and S. Weinberg, (World Scientific, Singapore, 1991), p. 269. 

[17] A. Macias, O. Obregon and J. Socorro, Int. J. Mod. Phys. A 8, 4291 (1993). 

[18] A. Csordas and R. Graham, Phys. Rev. Lett, 74, 4129 (1995). 

[19] O. Obregon, J. J. Rosales and V. I. Tkach, Phys. Rev. D 53, 1750 (1996). 

[20] P. D. D'Eath, Super symmetric quantum cosmology , (Cambridge University Press, Cambridge, 
England, 1996), p. 252. 

[21] H. Garcfa-Compean, O. Obregon and C. Ramirez, Phys. Rev. Lett. 88, 161301 (2002) 
[arXiv:hep-th/0107250j . 

[22] J. B. Hartle, in Quantum Cosmology and Baby Universes, eds. S. Coleman, J. B. Hartle, T. 

Piran and S. Weinberg, (World Scientific, Singapore, 1991), p. 65. 
[23] J. J. Halliwell, J. Perez-Mercader, and W. H. Zrek (eds.), The physical Origins of Time 

Asymmetry, (Cambridge University Press, 1994). 
[24] J.C. Thorwart, Quantum Cosmology and the Dechoherent Histories Approach to Quantum 

Theory, Ph.D. thesis, University of London, (2002). 
[25] J. J. Halliwell, Contemp. Phys. 46, 93 (2005). 

[26] J.J. Halliwell, Phys. Rev. D 36, 3626 (1987); Phys. Lett. B 196, 444 (1987). 
[27] T. P. Singh and T. Padmanabhan, Ann. Phys. (N.Y.) 196, 296 (1989); Class. Quant. Grav. 
5, 903 (1998). 

[28] H. Kodama, in Fifth Marcel Grossman Meeting, proceedings, Perth, Australia, 1988, edited 

by D.G. Blair and M.J. Buckingham (World Scientific, Singapore, 1989). 
[29] E. Calzetta and B.L. Hu, Phys. Rev. D 40,380 (1989). 

34 



[30] S. Habib and R. Laflamme, Phys. Rev. D 42, 4056 (1990). 

[31] M. Hillery, R.F. O'Connell, M.O. Scully and E. P. Wigner, Phys. Rep. 106, 121 (1984). 

[32] S. Habib, Phys. Rev. D 42, 2566 (1990). 

[33] R. Laflamme and J. Louko, Phys. Rev. D 43, 3317 (1991). 

[34] J. J. Halliwell, Phys. Rev. D 46, 1610 (1992). 

[35] H. Weyl, Group Theory and Quantum Mechanics, (Dover, New York, 1931). 

[36] E. P. Wigner, Phys. Rev. 40, 749 (1932). 

[37] A. Groenewold, Physica 12, 405-460 (1946). 

[38] J. E. Moyal, Proc. Camb. Phil. Soc. 45, 99 (1949). 

[39] J. Vey, Comment. Math. Helvet. 50, 421 (1964). 

[40] F. Bayen, M. Flato, C. Fronsdal, A. Lichnerowicz and D. Sternheimer, Ann. Phys. Ill, 61 

(1978); Ann. Phys. Ill, 111 (1978). 
[41] M. De Wilde and P.B.A. Lecomte, Lett. Math. Phys. 7 (1983) 487; H. Omori, Y. Maeda and 

A. Yoshioka, Adv. Math. 85, 224 (1991). 
[42] B. Fedosov, J. Diff. Geom. 40, 213 (1994); Deformation Quantization and Index Theory 

(Akademie Verlag, Berlin, 1996). 
[43] M. Kontsevich, "Deformation Quantization of Poisson Manifolds I". |q-alg/9709040|; Lett. 

Math. Phys. 48, 35 (1999). 
[44] D. L. Wiltshire, An introduction to quantum cosmology. [gr-qc/ 0101003|. 
[45] F. Antonsen, "Deformation quantisation of gravity," |arXiv:gr-qc/9712012|; "Deformation 

quantisation of constrained systems," |arXiv:gr-qc/9710021| . 
[46] C.K. Zachos, "Deformation Quantization: Quantum Mechanics Lives and Works in Phase 

Space", Int. J. Mod. Phys. A 17 (2002) 297. | hep-th/0110114| . 
[47] A.C. Hirshfeld and P. Henselder, Am. J. Phys. 70, 537 (2002). 

[48] G. Dito and D. Sternheimer, "Deformation Quantization: Genesis, Developments and Meta- 
morphoses", Deformation Quantization (Strasbourg 2001) Lect. Math. Theor. Phys. 1 Ed. de 
Gruyter, Berlin, IRMA (2002) pp. 9-54. 

[49] D. Sternheimer, "Quantizations and quantizations and introductory overview", Contemp. 
Math. 462, 41,(2008). 

[50] N. Seiberg and E. Witten, JHEP 9909, 032 (1999). |arXiv:hep-th/9908142| . 

[51] R.L. Stratonovich, Sou. Phys. JETP 31, 1012 (1956); A. Grossmann, Commun. Math. 

35 



Phys.48, 191 (1976); J. M. Gracia Bondi'a and J. C. Varilly, J. Phys. A: Math. Gen. 21, 
L879 (1988), Ann. Phys. 190, 107 (1989); J. F. Carihena, J. M. Gracia Bondi'a and J. C. Var- 
illy, J. Phys. A: Math. Gen. 23, 901 (1990); M. Gadella, M. A. Martin, L. M. Nieto and M. 
A. del Olmo, J. Math. Phys. 32, 1182 (1991); J. F. Plebahski, M. Przanowski and J. Tosiek, 
Acta Phys. Pol. B27 1961 (1996). 
[52] W. I. Tatarskii, Usp. Fiz. Nauk 139, 587 (1983). 

[53] H. Garci'a-Compean, J. F. Plebahski, M. Przanowski and F. J. Turrubiates, Int. J. Mod. Phys. 

A 16, 2533 (2001). [arXiv:hep-th/9909206| . 
[54] J. Luoko, Class. Quantum. Grav, 4, 581 (1987). 
[55] V.A. Rubakov, Phys. Lett. B 214, 503 (1982). 
[56] N. C. Dias and J. N. Prata, J. Math. Phys. 43, 4602 (2002). 
[57] N. C. Dias and J. N. Prata, Ann. Phys. 321, 495 (2006). 

[58] R. Kantowski and R. K. Sachs, J. Math. Phys. 7, 443 (1967), D. M. Solomons, P. Dunsby and 

G. Ellis, Class. Quant. Grav. 23, 6585 (2006). 
[59] C. Misner, in Magic Without Magic: John Archibald Wheeler (Freeman, San Francisco, 1972). 
[60] R. Szmytkowski and S. Bielski, J. Math. Anal. Appl. 365, 195 (2010). 

[61] A. Erdelyi et al., Tables of Integral Transform, Bateman Manuscript Project (McGarw-Hill, 
New York, 1954). 

[62] T. Curtright, D. Fairlie and C. Zachos, Phys. Rev. D, 58, 025002 (1998). 
[63] A. Erdelyi et al., Higher Transcendental Functions, Bateman Manuscript Project (McGraw- 
Hill, New York, 1953). 
[64] J. E. Lidsey, D. Wands and E. J. Copeland, Phys. Rept. 337, 343 (2000). 
[65] J. Maharana, Phys. Rev. D, 65, 103504 (2002). 

[66] E. C. Titchmarsh, Eigenfunction Expansions, Second Edition, (Oxford at Clarendon Press, 

Great Britain, 1962). 
[67] S. Fulling, Phys. Rev. D, 7, 2850 (1973). 

[68] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, 

New York, 1980). 
[69] S. W. Hawking, Phys. Rev. D 37, 904 (1988). 
[70] S. W. Hawking, Mod. Phys. Lett. A 5, 453 (1990). 
[71] A. Lyons and S. W. Hawking, Phys. Rev. D 44, 3802 (1991). 

36 



[72] I. Galaviz, H. Garcfa-Compean, M. Przanowski and F. J. Turrubiates, Annals Phys. 323, 267 
(2008) |arXiv:hep-th/0612245| . 



37 



